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The Hubbard model with local on-site repulsion is generally thought to possess a superconduct¬ 
ing ground-state for appropriate parameters, but the effects of more realistic long-range Coulomb 
interactions have not been studied extensively. We study the influence of these interactions on su¬ 
perconductivity by including nearest and next-nearest neighbor extended Hubbard interactions in 
addition to the usual on-site terms. Utilizing numerical exact diagonalization, we analyze the signa¬ 
tures of superconductivity in the ground states through the fidelity metric of quantum information 
theory. We find that nearest and next-nearest neighbor interactions have thresholds above which 
they destabilize superconductivity regardless of whether they are attractive or repulsive, seemingly 
due to competing charge fluctuations. 


The role of strong electronic correlations on uncon¬ 
ventional superconductivity (SC) remains an important 
open question [T] . The two-dimensional single-band Hub¬ 
bard model with an on-site repulsive Coulomb interac¬ 
tion has been well-studied as the simplest effective model 
which can capture the influence of correlations on the 
electronic properties. It is generally accepted that it 
possesses a superconducting ground state with d-wave 
symmetry upon carrier doping (changing the chemical 
potential); however, only a few studies can provide an ex¬ 
act solution, including a thermodynamic limit calculation 
limited to asymptotically small coupling Hi- In the 
intermediate to strong coupling regime, studies include 
numerically exact simulations HU, t - J model calcula¬ 
tions [ZlIH], and exact diagonalization on finite clusters 

[mi]. 

While the on-site Coulomb terms in the doped Hub¬ 
bard model are sufficient for producing a superconduct¬ 
ing ground state, one would like to understand the effects 
of more realistic, long-range Coulomb interactions. Of¬ 
ten, Hubbard models are used to approximate real mate¬ 
rials such as cuprates, so parameter choices are typically 
renormalizations of their bare values. As such, it makes 
sense to consider both repulsive and attractive extended 
interactions of various magnitudes [T2|. Naively, repul¬ 
sive (attractive) nearest neighbor interactions are ex¬ 
pected to suppress (enhance) SC, since d-wave pairs com¬ 
prise spins on nearest neighbor sites. Conversely, attrac¬ 
tive Tzext-nearest neighbor interactions favor diagonally 
adjacent spins but seemingly not SC’s nearest neighbor 
configuration, whereas repulsive interaction may enhance 
SC by suppressing diagonal adjacency. Calculations of 
repulsive nearest-neighbor interactions find that it indeed 
suppresses SC, but this does not mean that SC is desta¬ 
bilized. Although some contend that arbitrarily small 
interactions will destabilize SC m, several calculations 
show that SC remains stable until interactions reach some 
threshold UlHISj- On the other hand, attractive nearest- 


neighbor interactions are found in Hartree-Fock calcula¬ 
tions to enhance SC above some threshold when consid¬ 
ered alone m- However, one should also consider spin 
and charge correlations driven by these interactions, as 
they may enhance [20| or compete [2T] with SC. 

Here, we apply the fidelity metric from quantum infor¬ 
mation theory to explore the influence of extended Hub¬ 
bard interactions on SC. This method has been used to 
investigate classical as well as quantum critical behavior 
in various systems [22H32] and is a clear test of SC i- 
m- utilizing numerical exact diagonalization, we find 
that nearest and next-nearest neighbor interactions have 
thresholds above which they destabilize SC regardless of 
whether they are attractive or repulsive. 

Theory. - The two-dimensional extended Hubbard 
model in real space comprises kinetic and interaction 
terms, depicted schematically in Fig. The kinetic 
portion of the Hamiltonian is 

+ C^j„Cicr) + t ^ (c|^C;ct + , 

(ij)a ({il))a 

( 1 ) 

where {cia) creates (annihilates) a fermion at site i 
with spin a. () and (()) denote nearest and next-nearest 
neighbor sites, resp. Throughout, we set t = 1. The 
interacting portion of the Hamiltonian includes extended 
interactions up to next-nearest neighbors: 

Hint =U'^c\.^c\^CiiCi'^ + V ^ 

i (ij)<r<r' 

+ V ^ ( 2 ) 

We determine the ground state wavefunction for the 
model using exact diagonalization. The Hamiltonian ma¬ 
trix was diagonalized using the Parallel Arnold! PACK- 
age (PARPACK) [33] , which provides an efficient method 
to obtain a set of properly orthonormalized eigenvectors 
and is more numerically stable compared with standard 
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FIG. 1. (Color online) The operations of the kinetic and in¬ 
teraction Hamiltonian terms on spins: t and t' are nearest and 
next-nearest neighbor hopping, resp, U is the onsite interac¬ 
tion, and V and V' are next-nearest neighbor interactions. 
The lattice is a 16B Betts cluster with periodic boundary 
conditions, with atoms at the intersections. 


Lanczos techniques. Due to the exponential growth of the 
Hilbert space with cluster size, we work with a conserved 
electron number (canonical ensemble), fix the total mo¬ 
mentum (translational symmetry), and work in the total 
spin 0 sector. In a past fidelity study of the 16B Betts 
cluster HB, the Hubbard model with U W was found 
to have a superconducting ground state when doped with 
4 holes (25% filling) and with periodic boundary condi¬ 
tions (Fig. [^. Therefore, we study this cluster for the 
effects of extended interactions on SC. 

Having determined the ground state, one can under¬ 
stand the physical properties of the model using corre¬ 
lation functions. In our study, we are interested in the 
superconducting correlator, and the spin and charge cor¬ 
relations are also basic for our understanding. For d-wave 
SC, the pair density matrix 0 [HI [Ml isg between sites 
i,j at any separation distance in the ground state |4)o) is 
Pfp = where A = = 

CifCji + CjiCii- The charge density matrix is _ 

(^olPiPjl^o) - (^olPil^’o)^, where pi = The 

spin density matrix is where Si = 

Specifically, these correlators are used to construct 
charge and spin structure factors, which we analyze as 


a function of wave vector. 

fV(q) = (1/iV) ^ p^^harge^^q.r, ( 3 ) 

3 

5(q) = (l/iV)^PofV‘’--^ (4) 

i 

where N is the number of sites, q is the wave vector, and 
Yj is the position of site j. A structure factor to measure 
relative SC is similarly constructed with but it is 
preferable to use an absolute measure of the presence of 
SC. 

To investigate SC in Hubbard models, we calculate the 
fidelity, which was first used widely in quantum informa¬ 
tion to provide the criterion for distinguishing between 
quantum states [Ml ETj. This quantity has now been 
utilized increasingly in general physics for investigating 
classical as well as quantum critical behavior in vari¬ 
ous systems [2^32| . This metric measures the ground 
state wavefunction’s rate of change and shows whether 
the ground state adiabatically connects to a supercon¬ 
ducting state. Recently, it was used to demonstrate SC 
in the t-J 0(10] and on-site Hubbard models m- 

This technique involves adding an additional term to 
the Hamiltonian: 

Hsc = -^ZPlPp (5) 

The full Hamiltonian has the form H(Asc) = Ain + Ant + 
Ac- The fidelity metric is defined as 

g(Asc,<5A) = (2/fV)(l - P(Asc,<^A))/A^ (6) 

where the fidelity F = (4>o(Asc)l‘J’o(Asc+<5A)). For finite- 
size systems as in this work, a change in the ground 
state’s symmetry occurs via a quantum critical crossover 
(QCC) that appears as a broad peak in g. This provides a 
signal of strong SC: The ground state of H{Xsc) at large 
Age is expected to possess long-range SC with a correla¬ 
tion length at least of order the cluster size. If g displays 
no QCC as one reduces AgC) then this SC persists in the 
Hamiltonian of interest down to Age = 0. 

To substantiate the results of a fidelity analysis, we 
evaluate the ratio between the largest and next largest 
eigenvalues of the density matrices Pq^ and re¬ 

ferred to as R ratios [MIIIMIIM]- These eigenvalues are 
proportional to the structure factors for all values of q 
(Eq. [^and following text). Therefore, for charge correla¬ 
tions, if Pcharge > 1, there is a q for which the correlations 
are strongest, and Pcharge is the factor by which they are 
larger than the next strongest correlations. If Pcharge = 1, 
then there is no dominant q. For SC, only the familiar 
q = 0 correlations may be present in our model, so this is 
the dominant correlation for cases where Pgc > 1. Like 
correlators, the magnitude of P > 1 is not meaningful by 
itself, but its trends as a function of interaction parame¬ 
ters provide confirmation for the fidelity results. 
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FIG. 2. (Color online) The semi-logarithmic plot of fidelity 
metric vs. Asc for repulsive and attractive (a) nearest- 
neighbor and (b) next-nearest neighbor interactions on su¬ 
perconductivity with 4 doped holes. The insets show the fi¬ 
delity metric over a larger range of Asc to confirm that it is 
smooth without peaks for Asc > 0.5. Also shown are Rsc 
ratios for (c) nearest-neighbor and (d) next-nearest neighbor 
interactions and 7?charge density matrix eigenvalue ratios for 
the same (e,f). 


Results. - We show the effects of nearest and next- 
nearest neighbor interactions - both repulsive and attrac¬ 
tive - on superconductivity. The SC is clearest for the 
parameters U = 8,t' = -0.3 so we employ the same 
set here to analyze the effects of extended interactions on 
SC. To stabilize the numerics, steps of SX = 10“^ prove 
to be sufficient. 

Figure l^a) shows the fidelity metric for nearest neigh¬ 
bor interactions. For large AsC; the state is expected 
to to possess long-range SC with a correlation length at 
least of order the cluster size. Starting with a strong 
pairing field, Asc = 5, and decreasing AsCj the fidelity 
metric g increases monotonically (see inset of Fig. ©a)). 
This indicates an increased rate of change in the ground 
state as it approaches a potential QCC. For V = 0, g con¬ 
tinues to increase monotonically and does not peak down 
to Asc = 0; The ground state is still superconducting at 
Asc = 0, as found in the previous study El- 

Adding nearest neighbor interactions, we find that at¬ 
tractive and repulsive interactions have thresholds above 
which they suppress SC. The range of interactions shown 
is chosen to reach or cross these thresholds. For the re¬ 
pulsive case, y = 1 is similar to F = 0 but for V = 2, the 
QCC peak, which defines Ag^, is at 0.12; at Asc = 0 the 


system is no longer superconducting. Note that although 
this peak includes a sharp jump, the ground state’s total 
momentum remains at zero for all Asc shown. For at¬ 
tractive V < 0, g increases with decreasing Asc even for 
interactions as large as 1^ = -4 but with a much smaller 
slope than for smaller attractive interaction. For larger 
interactions Ag^ > 0 as SC is destabilized (not shown). 
This contradicts the intuition that attractive and repul¬ 
sive interactions have opposite effects on SC. 

Similarly for next-nearest neighbor interactions in Fig. 
Sb), we find that attractive and repulsive interactions 
have destabilization thresholds. For \V'\ = 1, g increases 
with decreasing Asc (see inset of Fig. ib)), indicating 
that the system remains superconducting. However, for 
larger V' = -2,2.2, one observes Ag^ > 0; the system 
is no longer superconducting, although the total momen¬ 
tum is unchanged. (At V = 2 the peak is still at Ag^ < 0.) 
Contrary to intuition, nearest and next-nearest neighbor 
interactions all destabilize SC beyond some threshold, re¬ 
gardless of whether they are attractive or repulsive. All 
of these destabilizing cases with AgQ > 0 are confirmed 
by level crossings in the ground state energy at Ag^ (not 
shown). It is possible that there are competing phenom¬ 
ena stabilized by these interactions. 

The fidelity results for nearest and next-nearest neigh¬ 
bor interactions are further confirmed by i?sc ratios in 
Figs. [^c,d), resp. In the cases where V or V desta¬ 
bilize SC and the fidelity metric peaks at Ag (2 > 0, 
Rsc drops sharply as Asc crosses below Xg^, corrobo¬ 
rating the suppression of SC. Additionally, the observa¬ 
tion that Rsc = 1 at Asc = 0 clearly indicates a non¬ 
superconducting system. Conversely, for the V, V' cases 
with stable SC, Rsc evolves smoothly for Asc > 0 just 
like the fidelity metric, implying no destabilization of SC 
in this range. 

To explore the possibility of phenomena competing 
with SC, we consider charge fluctuations (CF), which 
are analyzed via the i?charge ratios for nearest and next- 
neighbor interactions in Figs. [^e,f), resp. In all cases 
where SC is destabilized, as Rsc drops sharply at AgQ, 
.Rcharge increases sharply. This suggests that the increas¬ 
ing CF destabilizes SC. Consistently, in cases of stable 
SC, both Rsc and i?charge evolve smoothly. 

To find the wave vector of the CF for nearest neigh¬ 
bor interactions. Fig. [^a) displays the charge structure 
factor fV(q). At H = 2, where SC is destabilized, N{q) 
clearly peaks for CF with wave vector q = (tt, tt). For the 
other values of V that do not destabilize SC, (tt, tt) CF is 
suppressed. V <Q favors (f ,0) CF, which coexists with 
SC until it is destabilized for V < -A (not shown). Both 
attractive and repulsive interactions favor CF that seems 
to destabilize SC, but they have different wave vectors. 

We also consider the interplay between the spin struc¬ 
ture factor S'(q) in Fig. |^b) and superconductivity. At 
V = 0, the presence of SC is accompanied by (7r,7r) spin 
fluctuations (SF). These are suppressed along with SC 
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FIG. 3. (Color online) The (a) charge and (b) spin structure 
factors for all wave vectors q in the cluster at which they can 
be distinct in our model, for nearest neighbor interactions, 
and (c,d) the same for next-nearest neighbor interactions. 

at = 2 in favor of CF, whereas they increase for 1^ < 0 
where SC survives. It is CF, not SF, that appear to 
compete with SC at larger values of V. 

For next-nearest neighbor interactions, the CF is found 
via the charge structure factor in Fig. §c). Aty' = -2, 
where SC is destabilized, the structure factor A^(q) peaks 
for CF with wave vector q = (7r,7r), while for V' = 2.2 it 
peaks at q = (7r,0). For the other values of V that do 
not destabilize SC, no CF dominates. Both the attractive 
and repulsive interactions favor CF that destabilize SC, 
but with different wave vectors. Conversely, in Fig. [^d), 
the spin structure factor shows that q = (7r,7r) SF follow 
the behavior of SC and are thus suppressed at larger \V'\. 

Discussion and Conclusion. - Utilizing numerical ex¬ 
act diagonalization, we analyzed the signatures of super¬ 
conductivity through the fidelity metric, R ratios and 
structure factors. Consistent with intuition and past 
studies [I4lll8| . repulsive nearest neighbor interactions 
destabilize SC beyond some threshold seemingly due to 
competing CF. However, we find that this trend holds 
true for all nearest and next-nearest neighbor interac¬ 
tions, whether repulsive or attractive. 

These trends can be understood by simple considera¬ 
tions. The nearest neighbor repulsion and next-nearest 
neighbor attraction do not favor neighboring charge (see 
Fig. 0 and therefore drive (tt, tt) CF, consistent with 
lower dimensional na I3S] and weak coupling studies 
[131101 m]. The diagonally adjacent charge favored by 
(tt, tt) CF is incompatible with the nearest-neighbor spins 
of d-wave SC. The next-nearest neighbor repulsion favors 
(7r,0) CF for larger interactions. Even this wave vector 
may be incompatible with SC because it is somewhat 
restrictive on neighboring spins. Even if the obtained 


wave vectors may be subject to finite size limitations of 
the cluster, they are consistent with their driving inter¬ 
actions and other studies |42H45] . 

Even attractive nearest neighbor interactions, which 
may be expected to favor d-wave SC, will destabilize 
it for large enough interactions due to CF. This can 
be understood by considering the strong coupling limit, 
t/U, t/U ^ 0. For attractive \V\ » U, the electrons 
will gather as closely as possible. In our 16-site model 
with twelve electrons, the optimal configurations have six 
doubly-occupied sites with at most seven nearest neigh¬ 
bor bonds, by occupying a 2 x 3 block for example. The 
structure factors of these configurations peak at A^(|,0) 
(or equally N{0,^)) because their periodicity equals the 
size of the lattice. Their energy is A[y |»(7 = 28U + 6U. 
Conversely, in the \V\ « U limit, the sites are rarely 
doubly-occupied. Twelve out of sixteen singly-occupied 
sites have between sixteen and twenty occupied bonds, 
resulting in an energy E|y |«(7 between 16U and 20U. Set¬ 
ting A|yj »(7 = E\v\«Uj these energies cross at V between 
-U12 and -3f7/4, consistent with V = -C//2 at which 
we obtained the destabilization of SC by CF. Further¬ 
more, if such a strong coupling picture is appropriate in 
the thermodynamic limit, one may expect some thresh¬ 
old \V\ < U above which charge gathers together and SC 
is destabilized. 

In this work, we are able to access a broad phase space 
of extended interactions in a well-controlled manner via 
numerical simulations. The results are suggestive for 
the rich phase diagrams of strongly correlated materials, 
which may have varying effective extended interactions 
and screening lengths. These can give rise to or influence 
charge, spin, and superconducting behavior. 
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